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The renormalization group and operator product expansion are applied to the model of a passive 
scalar quantity advected by the Gaussian self-similar velocity field with finite, and not small, corre- 
lation time. The inertial-range energy spectrum of the velocity is chosen in the form E(k) oc fc 1_2e , 
and the correlation time at the wavenumber k scales as k~ 2+v . Inertial-range anomalous scaling for 
the structure functions and other correlation functions emerges as a consequence of the existence 
in the model of composite operators with negative scaling dimensions, identified with anomalous 
exponents. For 77 > e, these exponents are the same as in the rapid-change limit of the model; for 
77 < e, they are the same as in the limit of a time-independent (quenched) velocity field. For s = r\ 
(local turnover exponent), the anomalous exponents are nonuniversal through the dependence on a 
dimensionless parameter, the ratio of the velocity correlation time and the scalar turnover time. The 
universality reveals itself, however, only in the second order of the e expansion, and the exponents 
are derived to order 0(e 2 ), including anisotropic contributions. It is shown that, for moderate n, 
the order of the structure function, and d, the space dimensionality, finite correlation time enhances 
the intermittency in comparison with the both limits: the rapid-change and quenched ones. The 
situation changes when n and/or d become large enough: the correction to the rapid-change limit 
due to the finite correlation time is positive (that is, the anomalous scaling is suppressed), it is 
maximal for the quenched limit and monotonically decreases as the correlation time tends to zero. 

PACS number(s): 47.27.-i, 47.10. +g, 05.10.Cc 



I. INTRODUCTION 



In recent years, considerable progress has been achieved in the understanding of intermittency and anomalous 
scaling of fluid turbulence. The crucial role in these studies was played by a simple model of a passive scalar quantity 
advected by a random Gaussian field, white in time and self-similar in space, the so-called Kraichnan's rapid-change 
model [1|. There, for the first time the existence of anomalous scaling was established on the basis of a microscopic 
model [2\ and the corresponding anomalous exponents were calculated within controlled approximations |]-|(| and a 
systematic perturbation expansion in a formal small parameter |Q] . Detailed review of the recent theoretical research 
on the passive scalar problem and the bibliography can be found in Ref. ||. 

Within the approach developed in Refs. nontrivial anomalous exponents are related to "zero modes," that is, 

homogeneous solutions of the closed exact differential equations satisfied by the equal-time correlation functions. In 
this sense, the rapid-change model appears "exactly solvable." 

In a wider context, zero modes can be interpreted as statistical conservation laws of the particle dynamics ||. The 
concept of statistical conservation laws appears rather general, being also confirmed in numerical simulations by Refs. 
0,0, where the passive advection in the two-dimensional Navier-Stokes velocity field jl0| and a shell model of a 
passive scalar |nj were studied. This observation is rather intriguing because in those models no closed equations 
for equal-time quantities can be derived due to the fact that the advecting velocity has a finite correlation time (for 
a passive field advected by a velocity with given statistics, closed equations can be derived only for different-time 
correlation functions, and they involve infinite diagrammatic series). 

One may thus conclude that breaking the artificial assumption of the time decorrelation of the velocity field is the 
crucial point Jl(],|ll| . 

An important issue related to the effects of the finite correlation time is the universality of the anomalous exponents. 
It was argued that the exponents may depend on more details of the velocity statistics than only the exponents 77 and 



e [[12]. This idea was supported in Refs. 13 where the case of short but finite correlation time was considered for 
the special case of a local turnover exponent. In those studies, the anomalous exponents were derived to first order 
in small correlation time, with Kraichnan's rapid-change model Jll| or analogous shell model for a scalar field JTJ] 
taken as zeroth order approximations. The exponents obtained appear nonuniversal through the dependence on the 
correlation time. 
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In Ref. M and subsequent papers |15|-|19[], the field theoretic renormalization group (RG) and operator product 
expansion (OPE) were applied to the rapid-change model and its descendants. In that approach, anomalous scaling 
emerges as a consequence of the existence in the model of composite operators with negative scaling dimensions, 
identified with the anomalous exponents. This allows one to construct a systematic perturbation expansion for the 
anomalous exponents, analogous to the famous e expansion in the RG theory of critical behavior, and to calculate 
the exponents to the second [Q,[l5) and third |l6) orders. For passively advected vector fields, any calculation of the 
exponents for higher-order correlations calls for the RG techniques already in the lowest-order approximation |fnj| . 

Besides the calculational efficiency, an important advantage of the RG approach is its relative universality: it is 
not related to the aforementioned solvability of the rapid-change model and can also be applied to the case of finite 



correlation time or non-Gaussian advecting field. In Ref. |18| (see also |19[ for the case of compressible flow) the RG 



and OPE were applied to the problem of a passive scalar advected by a Gaussian self-similar velocity with finite (and 
not small) correlation time. The energy spectrum of the velocity in the inertial range has the form E(k) oc k l ~ 2e , 
while the correlation time at the wavenumber k scales as k~ 2+r) . It was shown that, depending on the values of 
the exponents e and 77, the model reveals various types of inertial-range scaling regimes with nontrivial anomalous 
exponents. For 77 > e, they coincide with the exponents of the rapid-change model and depend on the only parameter 
2e — r), while for e > 77 they coincide with the exponents of the opposite ("quenched" or "frozen") case and depend 
only on e. 

The most interesting case is 77 = e, when the exponents can be nonuniversal through the dependence on the 
correlation time (more precisely, on the ratio u of the velocity correlation time and the turnover time of the passive 
scalar). In the field theoretic language, the nonuniversality of the exponents in this regime is a consequence of the 
degeneracy of the corresponding fixed point of the RG equations. It agrees with the findings of Refs. [ p^|JT^ ] since 
the borderline rj = e, including the "Kolmogorov" point 77 = e = 4/3, corresponds to the case of a local turnover 
exponent. It is also interesting to note that the same relation r\ = e for the boundary between the time-decorrelated 
and quenched cases is encountered in a model of passive advection by a strongly anisotropic flow, studied in Refs. 
p0| . It was argued in ]2l]] that the same boundary will be observed with very general assumptions on the velocity 
statistics. Although the possibility of the nonuniversality of anomalous exponents for 77 = e was demonstrated by 
the rigorous RG analysis, the practical calculation by Ref. J^j has shown that they appear universal (independent of 
u) to the first order in e and 77: in the one-loop approximation, the anomalous dimensions of the relevant composite 
operators depend on a combination of the model parameters (couplings) that remains constant along the line of the 
fixed points. This fact is rather disappointing, because it means that in the one-loop approximation it is impossible to 
judge how the finite correlation time affects intermittency, in particular, whether the anomalous scaling is enhanced 
or suppressed in comparison with the rapid-change or quenched limits. 

In this paper, we present the anomalous exponents to order 0(e 2 ) (two-loop approximation) for the most interesting 
case rj = e, including the exponents of the anisotropic contributions, and study their dependence on u. [It is not 
necessary to separately consider the cases 77 > e (77 < e) , because the corresponding exponents are the same as for the 
rapid-change (quenched) velocity field and can be obtained from the case r\ = e i n the limits u — *■ 00 (u — > 0)]. 



In Sec. |II|, we describe our model and its interesting special cases. In Sec. ^Ij , we briefly recall the field theoretic 
formulation, the RG and OPE approach to the model, and the 0(e) result for the anomalous exponents The 
results of the two-loop calculation are presented and discussed in Sec. |y|: in S ec. IV A , we give the ano malou s 



exponents to order 0(e 2 ) and then discuss them separately for the isotropic (Sec. IV B| ) and anisotropic (Sec. IV C| ) 
contributions. 

The main conclusion of the paper can be formulated as follows: the qualitative effect of the finite correlation time 
on the anomalous scaling depends essentially on the correlation function considered, the value of u, and the space 
dimensionality d. For the low-order structure functions and in low dimensions (d = 2 or 3), the inclusion of finite 
correlation time enhances the intermittency in comparison with the both limits: the time-decorrelated (u — 00) and 
time-independent (u = 0) ones. Although the anomalous exponents have a well-defined limits for u — > 0, they show 
interesting irregularities in the vicinity of the quenched limit: a rapid falloff when u — increases from zero, with 
infinite slope for d — 2, with a pronounced minimum for u ~ 1. On the contrary, the behavior in the region of large 
u is smooth, like for the shell model studied in [jl4| . For higher-order structure functions and large d, the anomalous 
scaling is always weaker in comparison with the rapid-change limit and the corresponding (positive) correction is 
maximal for u = and monotonically decreases to zero as u tends to infinity. 

II. THE MODEL 

The advection of a passive scalar field 6(x) = 0(t,x) is described by the stochastic equation 

\7 t 6 = v d 2 6 + f, V t = d t +vA, (2.1) 
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where dt = djdt, di = d/dxi, vq is the molecular diffusivity coefficient, d 2 is the Laplace operator, v(x) = {vi(x)} is 
the divergence-free (owing to the incompressibility) velocity field, and / = f(x) is an artificial Gaussian random noise 
with zero mean and correlation function 



(f(x)f(x'))=C(t-t', r), 



(2.2) 



The form of the correlator is unessential; it is only important that the function C in Eq. ( |2.2|) decreases rapidly for 
r ^> L, where L is some integral scale. The noise maintains the steady state of the system and, if C depends on the 
vector r and not only its modulus r = r, is a source of large-scale anisotropy. In a more realistic formulation, the 
noise is replaced by an imposed constant gradient of the scalar field; see e.g. Refs. [pP Jl^ , p^ , 22 ] . 

In the real pro blem , the velocity field satisfies the Navier-Stokes equation. Following Refs. 
for v(a;) in Eq. (2T) a Gaussian distribution with zero mean and correlator 



14,h&hftE2l, we assume 



(Vi(x)vj(x')) 



W) 1 



Py(k) D v (t - t', k) exp[ik • (x - x')], 



where Py (k) = <5y — kikj/k 2 is the transverse projector and the function D v (t, k) will be chosen in the form 



D v (t-t',k) = ^- 



2u k d - 2+2£ 



exp[-uj(k)(t - t')}, uj(k) = u is a k 2 ^ 1 . 



(2.3) 



(2.4) 



Here D and Uq are positive amplitude factors and the positive exponents e and rj play the part of small expansion 
parameters in the RG theory; see Refs. |l^Jl^]. It is also convenient to introduce the "coupling constant" go = Do /is 2 
(expansion parameter in the ordinary perturbation theory). The infrared (IR) regularization is provided by the sharp 
cutoff in all m ome nt um intergals from below at k — m with m ~ 1/L. 

The model (2J3), {2A) contains two special cases that possess some interest on their own: in the limit uq — *■ oo, 
So = So/ u o = const we arrive at the rapid-change model, 

D v {u, k) -> g' Q v S(t - t') k~ d ~ 2 '+\ (2.5) 

while the limit uq — > 0, q'q = go/2uo = const corresponds to the case of a quenched (time- independent) velocity field, 

D v (tu,k)^g'y Q k- d+2 - 2£ . (2.6) 

The latter case has a close formal resemblance with the well-known models of random walks in random environment 
with long-range correlations; see Refs. [p3|,p4|. 



III. RENORMALIZATION GROUP AND OPERATOR PRODUCT EXPANSION 



The RG theory of the model (2.1 )— ( 2.4 ) is prese nt ed in Refs. [ |l 8[Jl S| ] in detail; below we briefly recall only the 
necessary information. The stochastic problem (2.1)-(f2~4|) can be cast as a field theory with action functional 



S{0, 8', v) = -v J D- 1 v/2 + 0'D f 9'/2 + 0' [-V t + v d 2 } 9, 



(3.1) 



where 9' is an auxiliary scalar field and Df and D v are correlators (2.2) and (2.3), respectively. In Eq. (3.1), all the 
required integr atio ns over x = (t, x) and summations over the vector indices are understood. 

The model ( |3.l| ) is logarithmic for e = rj = 0; the ultraviolet (UV) singularities have the form of poles in various 
linear combinations of e and rj in the correlation functions. They can be removed by the only counterterm of the 
form 9' ' d 2 9, which is equivalent to the following multiplicative renormalization of the parameters go, uq and vq in the 
action functional (JO): 



i/ = vZ v , 



= QU 2 ^Z, 



go = gn 



u = ujj? 1 Z u , 



(3.2) 



where g, u, and v are the renormalized counterparts of the bare parameters, fj, is the reference mass in the minimal sub- 
traction (MS) scheme, which we always use in practical calculations, and Zi = Zi(g, u; d; e, rf) are the renormalization 
constants satisfying the identities 



Zg — Z v 3 , 



Z u — Z v . 

Fixed points of the corresponding RG equations are found from the requirement that the (3 functions 



(3.3) 



3 



f3 g = V^g = g[-2e - r\ + 3j v ], @ u = V^u = u[-r) + j v ], 7„ = 2? M In Z v (3.4) 
vanish. Here T>u is the operation [id^ for fixed bare parameters and the relations between [3 functions and the 



anomalous di men sion 7„ result from the definitions and the relation (3.3). 

From Eq. (3.4) the exact relation (3 g /g — 3f3 u /u = 2(rj — e) follows which shows that the j3 functions cannot vanish 



simultaneously for finite values of their arguments, except for the case r\ = e. Therefore, to find the fixed points 
we must set either u = oo or u — and simultaneously rescale g so that the anomalous dimension 7„ remain finite. 



These two options correspond to the two limits (2.5) and (2.6), so that the rapid-change and quenched cases are fixed 



points of the general model. The analysis shows that the former is IR stable (and thus describes the inertial-range 
asymptotic behavior) for r\ > e while the latter is stable for r\ < e. 

Most interesting is the case ij = e, when the (3 functions become proportional, and the set (3 g — (3 U = reduces to 
a single equation. As a result, the corresponding fixed point is degenerate: rather than a point, one obtains a line of 
fixed points in the g-u plane. They can be labelled by the value of the parameter u, which has the meaning of the 
ratio of the velocity correlation time and the scalar turnover time. 

Existence of the IR stable fixed points implies certain scaling properties of various correlation functions at scales 
larger than the dissipative length ~ g Q 1 ' 3e , I n particular, for the equal-time structure functions 

5„(r) = ([^,x)-^i,x / )n, r = x-x' (3.5) 

one obtains 

5 n (r) = D n/2 r n ^^F n (mr) (3.6) 

(odd structure functions are nontrivial if the correlation function (vf) is nonzero or if a constant gradient of the 
scalar field is imposed) . In the presence of anisotropy the scaling functions F n (mr) can be decomposed in irreducible 
representations of the SO(d) group. In the simplest case of uniaxial anisotropy (which is sufficient to reveal all 
anomalous exponents) one can write 

F n (mr) = Pi{z)F nl {mr), z = (n • r)/r, (3.7) 

where Pi(z) is the Ith order Gegenbauer polynomial (Legendre polynomial for d = 3) and n is a unit vector that 
determines the distinguished direction. 

The leading behavior of the functions F n i for mr (inertial range) is found from the corresponding operator 
product expansion and has the form 

i^cx(mr) A "', (3.8) 

where the "anomalous exponent" A n i is nothing other than the critical dimension of the irreducible traceless Ith rank 
tensor composite operator built of n fields 9 and minimal possible number of derivatives ]I§| ]. For I < n such an 
operator has the form 

d n e ■ ■ ■ d tl 9 (d i ed l 9) p + ■ ■ ■ , n = l + 2p. (3.9) 

Here the dots stand for the appropriate subtractions involving the Kronecker d symbols, which ensure that the resulting 
expressions are traceless with respect to any given pair of indices, for example, diddjd — 5ijdkddk0/d. We also note 
that the numbers n and I necessarily have the same parity, that is, they can only be simultaneously even or odd. 

For the most interesting case of the degenerate fixed point, the dimensions A„; are calculated in the form of series 
in the only independent exponent e = rj, that is, 



A ni =]>> fc A« (3.10) 

fc=i 

In the lowest order one obtains [H: 



A (i) _ -n(n-2)(rf-l) + A t (d+l) 
nl ~ 2(d-l)(d + 2) [ ' 

with A; = l(d + 1 — 2). For k > 2, the coefficients A^ depend not only on d but also on the parameter u, the ratio of 
the velocity correlation time and the scalar turnover time, that labels fixed points in the g-u plane (see above). 
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The reader not interested in the details of practical calculation can skip the end of this section and pass to the 

(2) 

result for A 1/. Calculation of the higher-order coefficients in the e expansions for the rapid-change model is presented 
in Refs. [ [15|Jlq ] in detail. Analogous calculations for the finite correlated case are more difficult in two respects. First, 
there are more relevant Feynman diagrams in the same order of perturbation theory (for zero correlation time, many 
diagrams contain closed circuits of retarded propagators and vanish). Second, and more important distinction, is that 
the diagrams for the finite correlated case involve two different dispersion laws: u> oc k 2 for the scalar and uj oc k 2 ~ v 
for the velocity fields. As a result, the calculation, as well as expressions for the renormalization constants, become 
rather cumbersome already in the lowest (one- loop) approximation; see Refs. [|l8| , [L9|| . 

The latter difficulty can be circumvented as follows. Careful analysis shows that in the MS scheme all the needed 
anomalous dimensions, 7„ from (3.4) and j n i = In Z n i, in contrast with the respective renormalization constants 



Z v and Z n i , are independent of the exponents e and n in the two- loop approximation (for the one-loop approximation 
this is obvious from the explicit expressions; see JlJ,|l^|). ft is thus sufficient to calculate them for any specific choice 
of the exponents e and rj that guarantees UV finiteness of the diagrams. The most convenient choice is rj = and 
arbitrary e: all the diagrams remain finite, the exponents in the aforementioned dispersion laws become identical, and 
the practical calculations drastically simplify and become feasible. 

To avoid possible misunderstandings, it should be emphasized that such an independence is not guaranteed by 
the renormalizability of the model. The renormalizability in the analytic regularization only guarantees that the 
renormalization scheme can be chosen such that the correlation functions, along with the coefficients /3 and 7 in the 
RG equations, will be analytic at the origin in the space of two complex variables e and r\ |2j|. We used another 
scheme in which the functions j3 and 7 are independent of e and rj in the first two orders, which does not exclude 
nonanalytic dependence on these parameters in higher orders. We expect that in the three-loop approximation 
nonanalytic constructions like (e + rf)/(e + 2ij) will indeed appear in the anomalous dimensions, in particular, due to 
the n eces sity to take into account UV finite parts of the two-loop diagrams (with our choice of the sharp IR cutoff in 
Eq. (2.4), the one-loop diagrams have no UV finite parts; cf. ]I6| for the rapid-change case). 

Thus, we conclude that the knowledge of the renormalization constants for rj = is sufficient to obtain the anomalous 
dimensions, (3 functions, coordinates of the fixed points, and the critical dimensions of composite operators for arbitrary 
values of rj and e, including the most interesting case rj = e, which we always discuss from now on. 



IV. ANOMALOUS EXPONENTS IN THE TWO-LOOP APPROXIMATION 



A. General expressions 



We have performed the co mp lete two- loop calculation of the RG functions ( |3.4|) and the critical dimensions ( 3.10| ) 
of the composite operators (3.9) for arbitrary values of n, I, d, and u and obtained the following expression for the 
second coefficient in expansion (3.1C ) : 



A 



(2) 



* (d~l) 2 {d + 2) 2 (d + 4) 

+ 9C [n(d + n)(d - 1) - k(d + 1)]} j 

with A; = l(d + I — 2). Here and below we denote 
(u-l-l/u)(d+l) (d+1) 



^2(d + 4)A [n(n - 2)(d - 1 ) + A ; ] + (n - 2)\&B [n(n - 4)(d - 1) + 3A ; ] + 



A = 
B = 
C = 



2(d + 2)(l-u) 2(d + A)(l - u)u(l + u) : 
(d+1) 



-F, 



1 



3(1 -u) 2 (d + 4) 
I rQ„,2 



u + 1 



{u + lf 
I 



9(l-u) 2 l 4 2 U 



2(u + l) J (u + l) 
u[2d- l + u(d-2)] 



(u + iy 



2ud(d + 2) 
(l-d)(l-u) 

—F 3 { - 



Au(d+1) 



-F, 



(u + l) 
f 



F 2 



2(u+l) 
4(d+l) 



J(u,d), 



4 V 4 

[d + 1 + 2u(d - 2)] 



-F, 



(u + l) 2 



1 



F 2 



(u + iy 



u 2 (d + i) 1 

(d + 4) 3 \4j ' (u + f)(d + 4) i " \2(u+l)J (u + l) 2 {d + 4y* \{u+\) 2 



(4.1) 

(4.2a) 
(4.2b) 

)- 

(4.2c) 



We also have denoted Fk(x) = F(l, 1; d/2 + k;x) for the hypergeometric series 

ab a{a+l)b(b+l) V 
F(a,b;c;z) = l + -z+ ^-^ ■- 
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The values of Fk entering into (4.2) can be related by the recurrent relation (x— 1)^2 (x) = x(d+ 2)Fz{x)/{d+£) — 1, 
but the resulting expressions loo k mo re cumbersome and we shall keep both Fi and F3 in the formulas. 

The quantity J(u,d) in Eq. (4.2a) can only be expressed in the form of a single convergent integral, suitable for 
numerical calculation: 



J(u, d) 



r(d/2) 



0Fr((d-i)/2) J (« 

z(l-z 2 )V 2 (l-u-z 2 ) 



1 (1 _ z 2 )d /2 

dz 



[2(1+ U )-Z 2 ]V2 



arctan 



z 2 (l - z 2 )ln 



f 4uz 2 , 

z[2{l + u)- z 2 ] 1 ' 2 
(1 + «- z 2 ) 



1 



z(u — 1 + 2z 2 ) arcsinz- 



(4.3) 



where r(- • •) is the Eu ler Gamma function 



The quantities (4.2), and hence the dimensions 
,( 2 ) 



if), have finite limits for 



and 



0. 



In the first limit, 
for 

general I). The Q(l /u) correction to the rapid-change limit can be found from the following asymptotic expressions 
for the coefficients ([4.2]): 



coincides with the known result for the Kraichnan's rapid-change model (see for I = and 2 and |T 



A = 



B = 



C = 



(d+1) 



(I + 2/u) + 0(I/u 2 ), 



2(d+2) 

(l + 2/„> + 0(l/A 



12 \4 



(d+1) (I 
9(d + 4) 3 ^4 



^1f 2 (- 

6 2 I 4 



2(d+l) F A 
9(d + 4) 3 ^4 



(d-2) 
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(4.4a) 



(4.4b) 



+ 0(l/u 2 ). (4.4c) 



The opposite case, u = 0, corresponds to the quenched (time-independent) velocity field. This case was extensively 
studied in connection with the so-called "random-random walks" (random walks in r andom environments); see the 
review paper ]23| ] and references therein. Our results for the function (3 g from (3.4) and the corresponding fixed 
point are in agreement with the two-loop results qu oted in Ref. ]2j| for the model of random-random walks. To our 
knowledge, the dimensions of composite opera tors (3.9) have not been studied in t hat context, and below we give 
the asymptotic expressions for the coefficients (12) sufficient to find the dimensions ( |4.l|) up to order 0(u) near the 
quenched limit: 



A 



B 



(d + 1) (3d + 4) 



2d(d+ 2; 
(d+1) u(d+l 



+ u(d+l) 



(3d + 2) 



In 2 



(d+3) 



3(d + 2) 3 
(d+I)(d 2 +4) u 



2d(d-2) d(d + 2) (d + 4) (d + 4)(d + 6) 
I „ fl\ 4 



Fa 



2(d + 3) 
(d + 4) 2 



(d + 4) 3 V2 



9d(d + 2) 



9 ,(l-2d)F 2 2 



d(d + 2) 
4(d + l) 



(d + 4) 



4(d+l)(d + 2) 2(d 2 + 2d + 4) 



d(d- 2) 



(d + 2) 



(4.5a) 



(4.5b) 



(4.5c) 



up to corrections of order 0(u 2 ). 

It is worth noting that the 0(u) terms in Eqs. ( 4.5a ) and ( 4.5c ) contain poles in (d — 2) and thus diverge for d = 2. 
Analysis shows that, for d = 2, the leading correction to the u = result is not analytical in u and has the form 
it In it. Formally, the singularity at d = 2 is explained as follows. Some of the two-loop diagrams contain "energy 
denominators" of the form (k + q) 2 + 0(it), where k and q are two independent integration momenta. The numerators 
contain factors oc (P^ (k)^^) 2 stemming from transverse projectors in the propagators. These factors suppress the 
singularity at k = — q, occurring in the denominators for u = 0, and ensure the existence of the integrals over k and 
q, but the "collincar" divergence at k = q occurs if the 0(u) to the denominators is taken into account. Physically, 
this divergence can be related to a strong resonant interaction between the excitations of the passive scalar field with 
the opposite momenta k = — q of equal mod uli in two dimensions. We shall see below that this singularity remarkably 
affects the behavior of the dimensions (4T) for the values of d much larger than d — 2. 

Many studies have been devoted to the analysis of the inertial-range turbulence in the limit d — > 00 P,|26|-p8[| . Our 
model has no finite "upper critical dimension," above which anomalous scaling would vanish. Like in the rapid-change 

128 



case l 
revea 



and, probably, in the Navier-Stokes turbulence 27 



the anomalous scaling disappears at d = 00, but it 
s itself already in the 0(I/d) approximation. Along with the results jl| for the scalar rapid-change model, where 
the 0(l/d) expression for the anomalous exponents were derived for any e, this confirms the importance of the large-d 
expansion for the issue of anomalous scaling in fully developed turbulence. 
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Straightforward analysis of the expressions (4.2) shows that, for d — > oo, one has A = 0(d°) [it is important here 
that J{u,d) = 0(l/d)), B = 0(d°) and C = 0(d), namely, 



( U + 2)(3 U + 2) 

36(w+l) 2 v ' 



(4.6) 



It then follows that for large d, the dimension (4.1) behaves as 0(\/d 2 ) and is completely determined by the only 
contribution with the coefficient C. This gives 



A (2) _ (n-2)(n-l)(u + 2)(3u + 2) 3 
" ; ~ A{u + l) 2 d 2 +U(l/d ). 



(4.7) 



The general expressions ( |4.1| ) , fl4.2| ) are rather cumbersome, and in the subsequent sections we shall separately discuss 
isotropic contribution (even n, I = 0) and anisotropic ones (general n, I ^ 0). 



B. Isotropic sectors 



Expression (4.1) simplifies for the most important case of the isotropic sector (even n and 1 = 0): 

n(n — 2) 



A 



(2) 

7)0 



(d-l)(d + 2) 2 (d + 4) 



2(d + A)A + 6(n - A)B + 9{d + n)C 



(4.8) 



This gives = in agreement with the exact result A20 = fjs}| . This means that the second-order structure 
function is not anomalous. The formal proof is based on certain Schwinger equation, which has the meaning of the 
energy conservation law; it is almost identical to the analogous proof for the Kraichnan model, given in Ref . . 
For the simplest nontrivial case n — 4 one obtains 



A2) 
*40 



8(2A + 9C)/(d- l)(d + 2) 2 



(4.9) 



that is, the quantity B does not enter into the result. For n > 6, all the coefficients (4.2) contribute to the result. 
In Fig. [l], we show the behavior of the quantity 



\i0 



v (2)| 



/n 3 



(4.10) 



for n — 4, 6, 8 and 20 (from below to above) as a function of u for several values of d. We have subtracted the value 
of the dimension for the rapid-change case, such that the curves approach zero as u — > 00, and divid ed the difference 
by n 3 , such that the results for different n's have the same order of m agnit ude [the quantity (4.1) is a third-order 
polynomial in n]. It is worth noting that, since the leading coefficient ( 3.11 ) is independent of u, it drops from the 
difference A n o — A n o|«=oo of the exact dimensions, and in the leading order 0(e 2 ) the latter is proportional to the 
quantity ( n introduced above. 

As one can easily see from Fig. El the qualitative behavior of £„ depends essentially on the values of n and d. 
For moderate n and d (e.g., n — 4, 6 and 8 for d — 2 and 3), finite correlation time enhances the intermittency 
(anomalous dimensions become more negative) in comparison with the both limits: the rapid-change (u — 00) and 
quenched (u — 0) ones. While the 0(l/u) correction leads to a smooth decrease of £ n for increasing 1/u from zero (in 
agreement with the numerical simulation of Ref. p3 for a shell model near the rapid-change limit), the rapid falloff 
of C„ is observed for increasing u from the quenched limit u = 0. As a result of the competition between these two 
effects, ( n is not a monotonous function of u and has a pronounced minimum in the interval between and 1. 

We recall that the slopes of the functions £„ at u = are i nfini te in two dimen sions for all values of n due to the 
presence of poles (d—2) in the 0(u) terms in Eqs. ( 4.5a ) and ( 4.5c ); see Sec. IV A| . For d > 2, the slopes become finite 
but they still remain very steep for d = 3 and lead to a rapid falloff of £„, as Fig. 0(b) shows. This fact also suggests 
that the quenched case, in contrast with the rapid-change one, can hardly serve as a good zero-order approximation 
in studying more realistic models of passive advection by the velocity field with finite correlation time. 

If n or d is large enough, the minimum becomes less pronounced, the behavior of C„ becomes more regular (e.g., 
n = 20 for d — 2) and eventually ( n becomes a monotonically decreasing function of u (n = 20 for d — 3). For such 
values of n or d, the function £ n approaches zero at u — > 00 from above. In other words, the 0(l/u) correction to the 
rapid-change limit suppresses the i nter mittency, in contrast with the case of moderate n. 

In the limit of large d, from Eq. 4.7 one easily obtains 
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(4.11) 



(note that the 0(n 3 ) term in this appr oximation disappears and A„g becomes only quadratic in n). One can see that 
for all values of n, the difference ( 4.11 ) is positive, decreases monotonically when u grows, and approaches zero from 
above when u — > oo. 

In Ref. [PL3| , the 0(1 /u) correction to the rapid-change case was derived by the zero-mode techniques in the limit 
of large d and arbitrary (not small) values of e, for the case of a local turnover exponent (e = rj). Although the 
anomalous exponents were shown to be nonuniversal (dependent on u) , our resu lts d isagree with Ref. |1| in two 
respects. First, due to the universality (independence of u ) of the leading term (3.11), the ratio fl4.11|) is of order 
0(e 2 ) and not 0(e). Second, the 0(l/u) correction in (4.11) is positive for all n, while, according to |13[|, inclusion of 
the finite correlation time makes the anomalous exponents more negative in comparison with the rapid-change limit 
for all n and e. It is not clear whether this disagreement can be explained by some distinctions between our model 
( |2.3p , ( |2.4l) an d the velocity ensemble employed in Ejj. It is possible to show, however, that any modification of the 
function ( |2.4| ) consistent with the RG analysis performed in Ref. jL8| and section III above, leads to an universal 
(independent of u) expression for the leading term in A^V, so that the 0(1 /u) correction to the zero-correlated limit 
remains of order 0(e 2 ). 

The changeover from the behavior typical to low spatial dimensions to the behavior described by Eq. (4.11) also 
produces interesting patterns, as illustrated by Fig. |l| for d = 5 and 10. 



C. Anisotropic sectors 



Let us tu rn to the analysis of anisotropic contributions in the structure functions (3.6), (|3.7|), descri bed by the 
dimensions (3. 1C) with I 0. We recall that such contributions appear in the inertial-range expression (3.7) if the 
forcing (2.2) is chosen to be anisotropic, or a constant gradient of the scalar field is imposed. 



An important property of the first-order result (3.11) is that for any fixed n, the quantity A^V increases monoton- 



ically with I |15|. One can say that the exponents, associated with tensor composite operators (3.9), exhibit a kind 
of hierarchy related to the degree of anisotropy: the less is the rank I, the less is the dim ension and, consequently, 
the more important is the cor responding contribution to the iner tial- range expression (3.7). The leading terms in the 
even structure functions (J^jfy are given by the scalar operators (3J3) with / = 0, that is, they are the same as in the 



model with isotropic forcing (we recall that n and / should be simultaneously even or odd). 

This behavior is in agreement with the existing phcnomcnological ideas, according to which the anisotropy intro- 
duced at large scales by the forcing (boundary conditions, geometry of an obstacle etc.) dies out when the energy 
is transferred down to smaller scales owing to the cascade mechanism J2S|. The hierarchy of anisotropic contribu- 
tions appears rather universal, being also observed for a vector (magnetic) field, advected by the Kraichnan velocity 
ensemble |3(J , for the scalar advected by the two-dimensional Navier-Stokes velocity field |u| and for the turbulent 
velocity field itself Q . 

Nevertheless, the anisotropy survives in the inertial range and reveals itself in dimensionless ratios involving odd 
structure functions. 



7Zk — S2k+i/S. 



fc+l/2 



(4.12) 



For a number of models it was shown that the skewness factor TZi decreases down the scales but slower than predicted 
by phenomenological theories f^|, while the higher-order odd ratios (hyperskewness IZ2 etc.) increase, thus sign alling 



of persistent small-scale anisotropy p^ , p^| , p0| ,^ l |. Due to the aforementioned hierarchy of the dimensions ( 3.11 ), the 
leading terms in the odd structure functions ( p.5|) in our model are determined by the vector operators (3.9) with 
1 = 1, and it is easy to check the abo ve st atements from the explicit expression (3.11). 



Of course, the 0(e 2 ) contribution^ 4.1) cannot change all the above properties of the dimensions A n ;, determined 



by their leading O(e) contribution (3.11), as far as small values of e are concerned. However, since the dependence on 
u occurs only in the 0(e 2 ) contribution, it should be taken into account if one wishes to discuss how finite correlation 
time affects the hierarchy of the dimensions or the behavior of the dimensionless ratios. 

Consider the effects of the finite correlation time on the hierarchy of the anisotropic contributions. To this aim, 
consider the difference of the coefficients (|4.1|) for a fixed value of n and two neighboring subsequent values of I: 



A (2) _ A (2) _ 



2(21 + d) IC n (d, u) 
(d-l) 2 (d + 2) 2 (d + 4) 



K n (d, u) = {2(d + A)A + 9(n - 2) [2B - (d + 1)C) }. 



(4.13) 
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(we recall that for a fixed n, a ll possible values of I are either even or odd, so that the subsequent values of I differ by 
2). It is clear from Eq. fl4.13| ) that the sign and the dependence on u of the whole expression is determined by the 
behavior of the function K, n (d, u) . 

In Fig. ||, we plot the quantity K n (d, u) as a function of u for n = 2, 4, 6 and 20 (from above to below) for d = 2 
(a) and d = 3 (b). The function is always negative for all the cases studied and increases monotonically with u. This 
behavior persists in the limit of large d, as follows from the asymptotical expr ession (4.7). 

We thus conclude that the 0(e 2 ) contribution in the exact dimension (3.10) "tries to cope" with the hierarchy, set 
by the 0(e) term, for all values of n, I, d and u; this effect is at its strongest for u = and weakens monotonically as 
u increases from to oo. 

Now let us t urn to the dimensionless ratios TZf. in ( 4.12 ). From the discussion below Eq. ( [1.12 ) and asymptotic 
representa tion ( |3.6| ) one obtains the power-like inertial-range as ympt otic expression TZk oc (mr) lX2k+1 - 1 with A 2 fc+i.i 
from Eq. (3.10) [we recall that in our model A2 i o = 0; see Sec. IVBfl . Due to the it-independence of the first-order 
answer (3.11), the 0(e) contribution in the exponent A.2k+i l = e(d + 2 — 4k 2 )/2(d + 2) + 0(e 2 ) coincides with its 
analog for the Kraichnan model; see Refs. Q for k — 1 and jl9| for general k. It completely determines the qualitative 
behavior of the quantities TZk- for k = 1 one has A3,i > and the skewness factor TZ\ decreases with mr, while for 
k > 1 one has A2fc+i.i < and the higher-order ratios TZk increase for mr — » 0. 

In Fig. ||, we show the behavior of the second-order correction A^ +1 15 obtained from the general formula Eq. (4.1) 
and divided by (2k + l) 3 like for the even dimensions, as a function of u for k = 1, 2, 3 and 4 (from below to above) 
for d — 2 (a) and d — 3 (b). 

One can see that the effect of the 0(e 2 ) correction on the inertial-range behavior of the ratios TZk is different for 
different d, k and u. In two dimensions, the corrections are negative for all u and moderate k: the decay of the 
skewness factor TZi for mr — > appears even slower than indicated by the 0(e) approximation, while the growth of 
the ratios with k > 2 becomes faster. 

In three dimensions, the correction is negative for k — 1 so that the decay TZ\ for mr — > is also slower than in the 
0(e) approximation. For k = 2, the correction is negative for small u (so that the growth of the hyperskewness factor 
IZ2 for mr — > is faster than in the first-order approximation), but it changes its sign for some finite value of u and 
the growth of H2 slows down. For k > 3, the corrections are negative for all u and the growth of the corresponding 
higher-order ratios TZk appears slower than predicted by the 0(e) expression. One thus may conclude that for d = 3, 
with the exception of the k — 2 case, the effect of the second-order term is opposite to the tendency set by the 
first-order approximation. 

(2) 

For k large enough and any d, the behavior of the quantities A 2I J +1 1 becomes similar to that of the even dimensions 
(2) £ 1 

A 2k discussed in Sec. IV B: they decrease monotonically as u increases, comparatively fast for small u (due to the 



1 2fe,0 

singularity in the slope for d = 2; see Sec. IV A) and rather slow when u becomes large enough. This follows from 

behaves as 0(n 3 ) for n — > 00, while its 



the fact that the ^-independent contribution in the gene ral expression A r 

/-dependent contribution behaves only as 0(n); see Eq. ( |4.l| ). 

(2) 

We also note that for moderate k, the quantities A 2fc+1 1 show a nonmonotonous dependence on u in the re gion of 
small u and in this respect they also resemble the even dimensions; see Fig. [I] and the discussion in Se 



sec. 



IV B 



V. CONCLUSION 



We have applied the RG and OFE methods to a simple model of a passive scalar quantity advected by the synthetic 
Gaussian velocity field with a given self-similar covariance with finite correlation time. The structure functions of 
the scalar field exhibit inertial-range anomalous scaling behavior, as a consequence of the existence in the model of 
composite operators with negative scaling dimensions, identified with anomalous exponents. 

For the special case of a local turnover exponent, the anomalous exponents are nonuniversal through the dependence 
on a dimensionless parameter u that has the meaning of the ratio of the velocity correlation time and the scalar turnover 
time. The universality reveals itself only in the second order of the e expansion, and we have derived the exponents 
to order 0(e 2 ), including anisotropic contributions. 

It is shown that, for isotropic contributions, the qualitative effect of finite correlation time depends essentially on 
n, the order of the structure function, and the space dimensionality d. For moderate n and d, finite correlation time 
enhances the intermittency in comparison with the both limits: the rapid-change (u = 00) and quenched (u = 0) 
ones. The 0(e 2 ) term shows a highly nontrivial behavior in the vicinity of the quenched limit: a rapid falloff when 
u = increases from zero, with infinite derivative at u — for d — 2, with a pronounced minimum for u ~ 1. This 
irregularity shows that the time-independent advecting field can hardly be a reasonable approximation in studying 
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more realistic models of passive advection by the velocity field with finite correlation time. The behavior near the 
opposite limit, u — oo, is smooth in agreement with the existing simulation for a shell model 

The behavior changes remarkably when n and/or d become large enough: the correction to the limit u = co due 
to finite correlation time is positive for all u (that is, the anomalous scaling is suppressed in comparison with the 
rapid-change case), it is maximal for u = and monotonically decreases to zero as u tends to infinity. 

In the anisotropic sectors, the 0(e 2 ) terms diminish the hierarchy revealed by the first-order terms for all values of 
the parameters n, I and d; this effect is maximal at u — and decreases monotonically with l/u. 

The effect of the 0(e 2 ) corrections on the inertial-range behavior of the dimensionless ratios involving odd-order 
structure functions depends on d. For d — 2 and moderate k these corrections are negative; the decay of the skewness 
factor TZi for mr — > is slower while the growth of the higher-order ratios IZk with k > 2 is faster than indicated 
by the 0(e) approximation by Refs. P]l8|. For d = 3, the effect is, for most cases, opposite to the tendency set by 
the first-order approximation: the decay of the skewness factor slows down as well as the growth of the higher-order 
ratios. 

Our analysis has been confined with the region of small e, where the results obtained within the e expansion are 
internally consistent and undoubtedly reliable [we recall again that, although the leading terms of the anomalous 
exponents are of order O(e), the leading terms in which the effects of finite correlation time occur are of order 0(e 2 )]. 
We do not discuss here the serious issue of validity of the e expansions for finite e = 0(1). One can think that, in 
our model, the natural region of validity of the e expansion is restricted by the value e = 1/2, where the velocity field 
acquires negative critical dimension (along with all its powers) and new IR singularities, related to the well-known 
sweeping effects, occur in the diagrams; see the discussion in Ref. |Q. (It should be noted, however, that such 
singularities do not necessarily lead to a changeover in the inertial-range behavior, as shown in [Q for the special case 
of the structure function 5 2 for u = 0). On the other hand, e = 1/2 can be regarded as the upper bound of the range 
of validity of the model itself: the lack of Galilean covariance becomes a serious drawback of the synthetic Gaussian 
velocity ensemble when the sweeping effects become important. The next important step should be the analytical 
derivation of anomalous exponents of a passive scalar advected by the Galilean covariant velocity field; this work is 
now in progress. 
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